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[57] ABSTRACT 

Radiologists may protect their subjects against unnecessary 
exposure to X-rays by use of X-ray opaque material. The 
resulting image contains relatively unexposed regions which 
comprise the shadow cast by the X-ray opaque material. 
This invention discloses a method for the automatic deter- 
mination of the location of the boundary between signal and 
shadow regions. Many hypotheses as to the location of the 
boundary are generated and the final decision as to the 
correct hypothesis is contingent on the results of a number 
of high-level tests applied to each. 
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FIG. 3 
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METHOD OF RECOGNIZING AN 
IRRADIATION FIELD 

FIELD OF THE INVENTION ^ 

Hie present inveDtioa is in the field of digital radiography. 
The inventian more specifically relates to a method of 
recognising an inadiation field in a system wherein a 
radiation image is recorded on a photostimulable phosphor 
screen by means of an irradiation field stop (i.e. means for lo 
partially shielding against irradiation). 

DESCaEOPnON OF PRIOR ART 

Radiologists may protect their subjects against uimeces- 15 
sary exposure to X-rays by use of X-ray opaque (*collima- 
tion') material. The material is placed in the path of the 
X-ray beam so as to shield those areas of the patient which 
are not regarded as diagnosdcally important Besides reduc- 
ing patient dose this technique has the additional advantage 20 
of limiting the amount of scattered radiation in the image 
field of interest. The regions of the resulting image which 
comprise the shadow cast by the X-ray opaque material 
(shadow regions, cfr. example in FIG. 1) are exposed only 
by scattered radiation. 25 

In computer-assisted radiology, the X-ray image is digi- 
tised as a two dimensional array of numbers, whose mag- 
nitude is related to the intensity of X-rays arriving at the 
detector, via the patient The values in the array may be 
rescaled, prior to display or printing, in order to maximise ^ 
the visual contrast in the area of diagnostic interest Such 
rescaling depends, in part, on the intensity histogram of the 
X-ray image. To achieve this function automatically, best 
results are obtained if the 'shadow' cast by the coltimation 
material is excluded from consideration, since it cannot 
contain any useful information, but would otherwise domi- 
nate the lower end of the intensity histogram. 

The presence of the collimation shadow region can cause 
an additional problem in the display of the image on fihn. 
Therein, the shadow region is relatively bright, and if ^ 
displayed unmodified, may impair diagnosis of subtle 
lesions due to dazzle, especially if the unexposed region is 
relatively large. 

The U.S. Pat No. 4,952,807 by Adachi addresses the first 
problem, wherein a method is disclosed for adjusting the 
image processing conditions by selecting only those pixels 
which correspond to the object or to the background, i.e. the 
'signal* region, under the assumption that the collimation 
material does not paitiy cover the object, i.e. the shadow ^ 
region does not touch the object of the image, but rather its 
surrounding background. In many kinds of examinations 
however, parts of the patient will be shielded from exposure 
to avoid useless X-ray load, as is the case with lumbar spine 
and shoulder examinations. ^ 

The U.S. Pat No.' 4,804,842 firom Nakajima also 
addresses Che first problem, wherein a method is disclosed 
for removing the lower part of the histogram corresponding 
to the pixels in the shadow region in order to adjust the 
image processing conditions. go 

However, techniques for histogram modification alone, do 
not address the second problem, i.e. that of dazzle. Its 
solution demands an explicit detection of the shadow region, 
as distinct &om the signal region. The offending shadow 
region may then be excluded firom the display process. In 65 
addition, the first problem of rescaling for optimum visual 
contrast can then proceed on the signal region alone, without 
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the complicating influence of the shadow region on the 
intensity histogram. . 

Previous work in the field of automatically delineating the 
signal/shadow regions has relied on piecemeal attempts at 
the solution. 

Many of the proposed methods are based on collecting 
candidate signal/shadow border pixels (i.e. edge pixels at the 
interface between signal and shadow regions) by threshold- 
ing on first differences along predetermined lines that might 
cross the signal/shadow border. A typical example of this 
approach is proposed by Funahashi in the U.S. Pat No. 
4,970,393. The signal region is recogttized as a region 
enclosed by pixels where the first difference exceeds some 
specified threshold, eitiier positive or negative. 

A refinement of this technique is proposed by Takeo in the 
European patent 0 342 379, where multiple candidate pixels 
along a searching Ime are ranked, and only the highest rank 
candidates on each searching line are considered in assem- 
bling the signal/shadow boundary. In both methods the 
underlying assumption is that the signal/shadow boundary 
has strong contrast everywhere. In many cases however, 
contrast along parts of the signal/shadow boundary is nearly 
zero, especially if the collimation material is partially cov- 
ering very dense parts of the body. So if the recognized 
signal/shadow boundary is assembled &om high difference 
pixels many gaps will occur in general. 

In the European patent 0 285 174 by Funahashi a straight-^ 
forward application of the basic Hough transform technique 
is proposed to delineate the signal/shadow boundary, which 
is assumed to be polygonal. The Hough transform is applied 
to prospectiv e edge points obtai^ byHifierenuaaon. The 
uiiderlyihg assumptions are that the signal/shadow boundary 
is piecewise linear, and has high contrast with respect to the 
contrast of line features within the signal region. However, 
experimental evidence shows that edges within the signal 
region, e.g. bone edges, firequently have higher contrast than 
the signal/shack>w boundary. As a consequence points in 
Hough transform space that correspond to straight bone 
boundaries may have higher scores than some of the actual 
signal/shadow borders, in which case the subsequent delin- 
eation procedure will fail. 

A diffoent approach is disclosed in the European patent 
0 360 231 by Funahashi. The image is subdivided into 
mutually adjacent non-overlapping tiles. A statistical param- 
eter is computed indicating the edge content of the tile, 
either by differentiation in image space, or in frequency 
space. Tlie tiles are accordingly classified as belonging to the 
relevant image region or not. 

We ran elaborate experiments and concluded that single 
criteria based on edge contrast, or local statistical parameters 
at both sides of a prospective signal/shadow boundary, suc^ 
as mean value, min, max, variance or spatial fiequency 
coefSdents are inadequate in delineating the irradiated field 
with a high-level of confidence,. say >98 percent. We foimd 
that the signal/shadow border contrast will vanish in some 
cases; anatomical structures may be straight and have high 
contrast edges; the image signal (which represents local 
dose) in the shadow region will occasionaly be higher than 
the signal in the unshielded (i.e. diagnostic) image region — 
in cases where the scattered radiation is high and some parts 
of the body are difficult to penetrate by radiation. Also we 
found that there is no essential distinction between image 
areas corresponding to very poorly penetrated bone, and 
areas corresponding to shielded parts. Both kinds of areas 
will be characterised by an extremely low mean value (i.e. 
dose), similar noise level (both quanmm limited), and by the 
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absence of other features (since the signal to noise ratio is 
very low). . 

In order to obtain a high success rate, a preferred method 
should consider many hypotheses as to the location of the 
signal/shadow boundary, and keep or reject each prospective 
signal-shadow boundary "as a whole", instead of keeping or 
rejecting individual components of a hypothesis, such as 



Multiple hypotheses as to the location of the signal/ 
shadow boundary are considered in the method disclosed by 
Tkiaka in the U.S. Pat No. 4,952,805. For each of the 
corresponding candidate signal regions, two signal histo- 
grams are computed: one for the pixels inside the candidate 
irradiation field and one for the pixels outside. The caxididate 
boundary which is characterised by the largest degree of 
interclass separation as determined from both histograms is 
selected to be the recognised irradiation field. The degree of 
interclass separation is defined in such a way that a high 
value will result if there is very little overlap between both 
Mstogranis. Hence this criterion will be effective only if the 
range of signal values within irradiation field are well 
separated &om the signal values in the shielded region, 
which docs not hold in many practical cases as we verified 
on the basis of digital radiographs from various examina- 
tions. 

Also, there is a practical limitation regarding the usability 
of the above mentioned approach, in that the number of 
hypotheses as to the locadon of the irradiation field may be 
prohibitively large. 

Previoiis work in general has focussed on methods in 
which finding the boundary of an irradiation field essentially 
relies on a very small number of mostly local criteria. 
Reviewing a wide variety of examinations confirms that 
individiial criteria will only wort in a very limited number 
of cases. A high success rate can only be achieved if many 
sources of evidence are combined, ranging from local image 
features to glc^al constraints. 

The behaviour of the radiologist, in placing the collima- 
don materia] aroimd the patient, may broadly be character- 40 
ised by a set of simple rules. This invention addresses only 
the use of collimation material which casts a shadow whose 
outline is piece-wise linear. More specifically, the boundary 
between signal, and shadow regions is assumed to be a 
rectangle of arbitrary location and orientatioa It is assumed 
that the signal region is on the interior, and that any part of 
the rectangular boundary may, or may not, have a projection 
onto the (rectangular) detector, and therefore the iinage 
array; see FIG. 1. 
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OBJECTS OF THE INVENTION 

It is an object of the present invention to provide a method 
of recognising an irradiation field in a system wherein a 
radiation image has been reconled by using an irradiation 
field stop (collimation material). 

It is a further object to provide a method for determination 
of the location of the signal/shadow boundary of an X-ray 
image, wherein the x-ray image is represented by a digital 
signal representation and wherein a signal/shadow boundary 
is the locus of interface between the signal and the shadow 
regions resulting from the obscuring by the collimation 
material. 

It is a still further object to provide such a method in a 
system wherein a radiation image is recorded in a photo- 
stimulable phosphor screen and wherein said radiation 
image, is read by scarming the screen with stimiilating 
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radiation and by detecting the light emitted upon stimula- . 
tion. 

It is still a ftuther object of the present invention to aid the 
process of contrast-enhancement of an X-ray image, by 
permitting such to perform its calculations based solely on 
the diagnostically useful signal part of the X-ray iinage, and 
to reduce the brightness in shadow region to avoid dazzle 
during diagnosis. 

STATEMENT OF THE INVENTION 

The present invention provides a method for the deter- 
mination of the location of the signal/shadow boundary in an 
X-ray image represented in a digital signal representation 
comprising the steps of: 

i) Extracting low-level primitives from the X-ray image 
{X(io)}. 

ii) Forming a reduced number of intermediate-level primi- 
tives from the low-level primitives, 

iii) Building hypotheses as to the location of the signal- 
shadow boundary from combinations of intermediate- 
level primitives, during which each combination is 
subject to intermediate-level tests, whereupon partial or 
complete hypotheses are rejected or accepted, 

iv) Performing high-level verification tests on each pro- 
spective signal-shadow boundary whereupon hypoth- 
eses arc rejected, or accepted at some cost, 

v) Selecting the prospective signal-shadow bouridary with 
die least cost 

The method of the present invention is a major departure 
from earlier methods in that, in general, it generates many 
hypotheses as to the location of the boundary. The final 
decision as to the correct prospective signal-shadow bound- 
ary is contingent on the results of a number of high-level 
tests applied to each. 

The method permits codification of expert-system type 
knowledge as to the geometry, topology, and statistics of the 
signal/shadow boundary, and associ^ed regions. 

The invention is, for example, applicable to a digital 
image signal obtained by exposing a photostimulable phos- 
phor screen to an image originating from a partially shielded 
exposure of an object to X-rays, by scanning the screen with 
stimulating irradiation, detecting the light emitted upon 
stimulation and converting the detected light into a digital 
signal representation. 

Commonly the digital signal is subjected to image pro- 
cessing such as contrast enhancement and signal-to-density 
m^jping before being applied to an output device for 
reproduction or display. The result of the process of contrast- 
enhancement of an X-ray image is enhanced by permitting 
such to perform its calculations based solely on the diag- 
nostically useful signal part of the X-ray image, and to 
reduce the brightness in shadow region to avoid dazzle 
during diagnosis. 

In the statement of the invention and in the description 
hereinafter the term **X-ray image" is used in the meaning of 
a two-dimensional digital image repesentation as an array of 
immbcrs, the magnitudes of which are directly related to the 
intensity of the X-rays arriving at a 2-dimensional physical 
detector such as a photo-stimulable phosphor screen. 

By "primitives" is meant features, such as edges, which j 
are extracted from die image as the result of low-level image I 
processing, cfr. FIG. 2. 

An "intermediate-level primitive" is a combination of 
(basic) primitives which may be considered collectively 
when bitilding a prospective signal-shadow boundary. Inter- 
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mediate-level primitives are generally fewer, and (spatially) 
larger than ate (basic) primitives. Each primitive is con- 
structed from a mutually exclusive set of pixels. Each 
intermediate-level primitive is constructed from a mutually 
exclusive set of (basic) primitives. 

By the term "signal region" is meant the spatial area of an 
image which is not obscured by collimation material firom a 
direct line-of-sight of the X-ray source, see FIG. lb. 

Hie "shadow region" is the area of the image which is 
obscured by the collimation material from a direct line-of- 
sight of the X-ray source; see also FIG. lb. 

The "signal-shadow boundary" is the locus of the inter- 
face between the signal and the shadow regions. Kg, the 
pixels at the interface between the signal sad the shadow 
regions. 

A low-level primitive can be a line— being a data stmc- 
ture to represent an edge, and associated data, including its 
parameterisation — or a tile — being a small rectangular area 
of pixel resulting ficom a tessellation of the image array. 

Steps i) and ii) of the method of the invention achieve a 
reduced set of primitives which can be assembled into 
hypotheses as to the location of the signal-region. 

A prospective signal-shadow boundary is defined as a 
possible solution posited by the method of the present 
invention, said solution being a segmentation (division into 
similar regions) of the original image into signal and shadow 
regions, e.g. a two-dimensional array of numbers isomor- 
phic with the X-ray image, with pixel labels for signal and 
shadow regions. 

Hypotheses are generated ftom the reduced (intermediate- 
level) set because, in general, the number of low-level 
primitives is too num^us to perform an exhaustive evalu- 
ation (step iv) of every possible combination (step iii). 
Specifically, if the low-level primitives are hnes. then tiie 
\ intermediate-level primitives would be line-clusters, die 
\ result, for instance, of a linear agglomerative clustering 
j process. 

I A line cluster is then one or more lines which are deemed 
' to be in some manner similar, and have been accordingly 
grouped together. 

E.g. one or more lines which have been grouped together 40 
by a linear agglomerative clustering process. In the example 
of FIG. 3, the lines lb and Ic may be clustercd together to 
form a line cluster {lb, Ic}. Similariy, in the same example, 
V/the likely dusters are {la}, {Id}, {Ic. 1/ 1^}. 

Alternatively, a low-level primitive, can be a tile, i.e. a 45 
small rectangular array of pixels resulting from a tessellation 
of the image array. In this case an intermediate-level primi- 
tive would be a tile-duster, being a group of spatially 
connected tiles having in conmaon a feature or an attribute. 
Step iii) combines the intermediate-level primitives into 

complete hypotheses posited by the system. 
The definidon of the term prospecdve signal-shadow 
boundary was given hereinbefore, by which term Is meant a 
possible solution regarding the shape and position of signal 
and shadow-regions. 

Such a solution is a segmentation or division of the image 
into signal and shadow regions. 

One of the possible representations of such a segmenta- 
tion is a two dimensional array of numbers, isomoq^ with 
the X-ray image, the values of which represent a unique 60 
label which identifies the corresponding region in the X-ray 
image as either signal or shadow region. 

A labd generally is a pixd-value which associates the 
pixel, i.e. the location in the image array, with a particular 
class (having a particular quality or attribute). 

In Older to limit the number of hypodieses to be submitted 
for evaluation (step iv). the method includes provision for 
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. tests which may be applied during the construction of the 
prospective signal-shadow boundary, which serve to detect 
and reject, at an early stage, an incorrect prospective signal- 
shadow boundary. 

Specifically, if the intermediate-level primitives are line- 
clusters for instance, then a possible test is on the mutual 
orthogonality of the line-clusters as putative sides of the 
signal/shadow border. Thus combinations of line-clusters 
which are not going to form a rectangle can be rejected at an 
early stege. 

Step iv) permits the evaluation of hypotheses successfully 
assembled from the intermediate-level primitives. A 
major feature of this invention is the generation of 
multiple hypotheses, where in each prospective signal- 
shadow boundary is a candidate for the complete 
signal/shadow boundary. 

The high-level tests permit evaluation of each prospective 
signal-shadow boundary, specifically of those attributes 
which can only be determined once a complete signal/ 
shadow boundary is prospective signal-shadow boundaryed. 
By such high-level tests, hypotheses at this level may be 
ramked (or possibly rejected) in order that a single candidate 
may eventually be chosen. 

Spedfically, the ranking of hypotiieses may be based on 
a scalar quantity representative of either the net positive or 
negative evidence in favour of, or against that prospective 
signal-shadow boundary. The distinction is immaterial: the 
highest ranking prospective signal-shadow boundary is that 
most like the idealised model. In the former case this will 
have the largest (positive score), in the latter: the smallest 
negative score. 

The evidence may be based on any of a number or 
combination of primitives, and less local features of the 
original X-ray image. For instance, negative evidence may 
be accumulated against a candidate signal/shadow boundary 
for each comer of the signal/shadow boundary according to 
(a monotonic function of) the magnitude of the difference of 
the angle at the comer from 90 degrees. Such a test would 
attempt to enforce the rectangular geometry of the modd. 
Broadly, tests may seek to enforce the geometry, topology, 
region statistics, and edge statistics required by the modd. 

Mthin these categories, other possible sources of evi- 
dence may involve tests for 
Edge statistics: 

Linearity: The straightness for each of the candidate 
signal/shadow boundary sides (however they were 
achieved). 

Existence: The accumulated evidence along the length of 
each side, computed &om an edge-image, such as that 
furnished by Carmy edge-^tection for instance. 

Continuity: The length of each candidate signal/shadow 

. boundary side fw which there exists no evidence from 
the edge-detection process. 

Explanation: The ratio of used evidence to the total 
evidence. In this context, the evidence is the number of 
edge-pixels. The used evidence is the immber\of edge 
pixels in the signal/shadow boundary. The unused 
evidence is the number of edge pixels which lie on 
edges outside of, but colinear with, those edges deem- 
ded to be in the candidate signal/shadow boundary. 
Geometry: 

Proportion: The aspect ratio of the rectangle (peiiiaps 
favouring rectangles with an aspect ratio similar to that 
of the X-ray image). 

Area: Hie area of the candidate signal-region (perhaps 
penalising very small rectangles). 
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Position: The placement of the centre of the candidate 
rectangle (perhaps favouring rectangles with their cen- 
tre near the centre of the X-ray image). 

Orientation: The orientation of the candidate rectangle 
within the image (perhaps favouring rectangles parallel 
with the edge of the X-ray image). 
Region statistics: 

Exposure: The mean, maximum, minimum (and the like) 
statistical measures derived from the (candidate) signal 
and shadow regions of the X-ray image (perhaps 
favouring low-exposure grey-levels in die shadow, 
versus signal regions). 

Homogeneity: The variance (and the like) statistics 
derived from the (candidate) signal and shadow regions 
of the X-ray image (perhaps favouring low-variance 
grey-levels in the shadow, versus signal regions). 
Tbpology: 

\^8ibility: The visibility of the candidate signal/shadow 
boundary, i.e. the degree to which the candidate signal/ 20 
shadow boundary is wholly contained within the image 
(perh^s favouring complete, rather than clipped 
boundaries — eg as in FIG. Id rather than FIG. Ic). 

Completeness: Hie number of visible rectangle sides. 
Hypodiieses consisting of just two *sidcs* (which are 25 
approximately orthogonal or parallel), can be consid- 
ered as partially visible rectangles. The scoring may, 
for instance, favour hypotheses with more, over fewer. 



10 



15 



The scalar result of each test may provide for the rejection 30 
of a prospective signal-shadow boundary by assuming an 
^;jpropriately small positive or large negative value, effec- 
tively denying that prospective signal-shadow boundary any 
chance of highest ranking. 

Step v) is just the selection of that prospective signal- 
shadow boundary which is ranked the highest of those 
surviving the higji-level tests of step iv). 

The method of the present invention can be applied for 
recognising an irradiation field in cases where at the expo- 
sure step an irradiation field stop has been used. 

It will be clear that the invention is not liinited to this 
particular application. It can for example also be used when 
the borders of irradiation fields have to be determined in an 
application wherein at least two images are recorded in a 
juxtaposed relation, a mode often referred to as *split 
screen'. 
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BRIEF DESCRIPTION OF DRAWINGS 

The method of the present invention as well as particular 
embodiments hereof will be explained hereinafter with 
reference to the accompanying drawings wherein 

FIG. 1 shows examples of location, orientation, and 
visibility of signal/shadow boundary, 

FIG. la shows an image with no collimation shadow 
regions, 

FIG. lb shows an image wherein two sides of the rect- 
angle are visible, (1) is the signal region, (2) are the shadow 
regions, 

HG. Ic shows an image with 3 sides of the rectangle 
visible, 

FIG. Id shows an image with all 4 sides of the rectangle 
visible, 

FIG. le shows an image wherein 4 sides are partially 
visible. 
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FIG. 1/ shows an image with 1 side visible, 

HG. 2a shows the original image, 

FIG. 2b shows the edge image. 

FIG. 3 is an example for line clustering. Lines 1^, and 1^ 
may be clustered. Similarly, lines 1^, Ip and 1^ may be 
clustered. 

FIGS. 4a-b are an example of edge tMrming, 
FIG. 5 are examples of vertex classifications for the 
central pixel in each 3x3 array, - 
FIG. 5a; isolated point 
FIG. Sb: terminal, 
FIG. 5c; edge-point, 
FIG. Sd: T-junction, 
FIG. 5e; X-junction. 

FIG. 6 is an example of linear extrapolatioiL Points and 
P^ arc the Hmiu of line l^. 

FIGS. 7c-c arc examples of closest zqiproach of two lines, 

FIGS. Sa~b are an example showing overlapping, and 
non-overiapping lines. 

HG. 9a shows an edge image, 

HG. 9b shows vertex image edge-pouits, 

HG. 9c shows projection of line-clusters. 

HG. lOa-h illustrates all legal combinations of line 
clusters. 

HG. 11 is a possible search tree for example of HGS. 9 
and 10, the letters denote the labelling in HG. 10, 

HG. 12 is an exposure and read-out apparatus. 

The following is a detailed description of each of the 
individual method steps of the present invention. 

Step i) Extraction of low-level primitives 

Tlie preferred low-level primitives are lines. They are 
derived from the X-iay image in two stages: edge-detection 
and segmentation. 
1.1 Edge Detection 

The purpose of this stage is to produce an edge-image: see 
for example HG. 9a. Our system incorporates a Canny 
edge-detection process. 

The following are definitions of some terms used in 
connection with a Canny edge detection process: 

Carmy edge-detection: A particular process for edge- 
detection comprising optional smoothing, followed by pro- 
duction of an edge-image wherein a pixel is labeUed as an 
edge-pixel if and only if the Sobel squared magnitude at that 
point is a local maximum along the (one-dimensional) 
trajectory which is normal to the Sobel oxientation. 

A Canny edge-image {E (r)} is formed by non-maximal 
suppression of the Sobel squared magnitude in the direction 
of the orientation of the image gradient, i.e. normal to die 
Sobel orientation: 



55 



60 



where 



1 if G(£)>G(r+^) a C(r- 
Ootheiwise 



65 



(l,-l)if«C(r) + it/4l -wOXnffl 

(l.l)ifie(i:)+nAH<ic« 

(0,l)ifli2(r)I^7c/8 



Let a mean a is bold and vectors r. represent the 
two-dimensional coordinates on the image plane (r=<ij) 
say). For an illustration, see HG, 2. 
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Smoothed image: An intrinsic image derived from the 
X-ray image wherein the high spatial frequencies are sup- 
pressed relative to the low spatial frequencies. E.g. the result 
of convolution of the X-ray image with a (two-dimensional) 
Gaussian. 

Sobel squared magnitude: The pixel-value G(i j) repre- 
senting the square of the gradient-magnitude of the image 
{X(iJ)} as determined by the following fomiula: 

Sobel orientation: The pixel-value Q (ij) representing the 
orientation of the tangent to the iso-gtadient-magnitude of 
the image {X(i j)}, as determined by the following formula: 

where atan2 (y^) is the arctangent C language function 
defined in [-71, n], 

Th& edge-detection process incoiporated in our system 
comprises four stages: 

1.1.1 Subsampling of the X-ray image 

The subsampling reduces the number of pixels of an 
image. For example: an original image obtained by reading 
an irradiation image stored in a photostimulable phosphor 
screen (cfir. FIG. 12) corr^rising approximately 5 million 
pixels, is reduced to an image of about 20000 12-bit pixels, 
preserving the aspect ratio of the original image. 

1.1.2 Canny Edge-Detection 

The subsampled image is passed to a Canny edge-detector 
to produce images: Sobel-squared magnitude {G(iJ)}, 
Sobel-orientation {Q(iJ)}, and edge-image {E(ij)}. The 
Gaussian smoothing is via a symmetric 5x5 convolution 
mask with weights calculated to effect a standard deviation 
of 035 pixels. 

1.1.3 Elimination of Weak Edges 

Edge-points with a low gradient magnitude indicate a 
weak edge. Studies have shown that the signal/shadow 
boundary is generally of high contrast, giving a large gra- 
dient magnimde. Very weak edges are therefore unlikely to 
form part of the signal/shadow boundary. Those with a 
Sobel-squared-magnitude below a threshold of 2000 are 
renwved from {E(ij)}, and correspondingly from {Q(ij)}, 
and {G(ij)}. 

1.1.4 Angle of Locus 

The orientation {0(i j)}, of the locus of edge-points in the 
edge-image {E(io)}, is calculated in three stages: edge- 
tMmiing, connectivity analysis, linear regressioa 

1.1.4.1 Thinning 

The edge-image is destructively thinned. The motivation 
is threefold 

1. The parameters of a regressed straight line, including 
the goodness-of-fit, are more accurately determined 
from the spatial distribution of edge-points. 

2. The various types of junction (TERMINAL. T-JUNC- 
TION, . . . ) are more easily detected. This is required 

_for comer-detection (below). . 
3- There then exists a close OKrespondence betwe^ the 
length of a straight line, and the number of (connected) 
edge-pixels of which it is made up. This is required for 
line-statistics (above). 
Thinning is achieved by removing edge-pixels which do 
not affect the local (3x3) coimectivity of other edge-pixels. 
See FIG; 4 for an example. 

1.1.4.2 Connectivity Analysis 

Connected edges are uniquely labelled by recursive edge- 
walking to create a label image {L(ij)}. All edge-pixels 
mutually cotmected share a common, but otherwise unique, 
label. 
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1.1.4,3 Linear regression 
In this connection the following definitions are first given: 
linear regression: The process by which an edge can be 
optimally parameterised according to the spatial loca- 
tion of its constituent edge-pixels. E.g. the determina- 
tion of the parameters of a line which are the radius to 
the normal (from the origin), and its orientation relative 
to the X-axis, by a least squares (of perpendicular 
distances) fit. 

Weighted linear regression: The process by which an edge 
can be optimally parameterised according to both the 
spatial location of its constituent edge-pixels, and the 
degree of certainty as to their correct classification 
Oabelling). E.g. the determination of the parameters of 
a line which are the radius to the normal (firom 
the origin), and its orientation Q^,^ relative to the 
X-axis, by ihc minimisation of a weighted least squares 
(of perpendicular distances) fit, the weighting being the 
Sobel squared magnitude. 

Let Pj^ be ihc positions of the edge-pixels (indexed by k) 
which comprise the line: 

and let the positional covariance matrix be F: 



given by 

for a,b in {x,y} 
then 



and 
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The eigenvalues of F are the sum of square errors 
(i.e. proportional to the variances) in each of the principal 
directions le along, and normal to, the best fit straight line: 

2Xt=rr(n±(rf(r)*-4/)tf;(r))^ 

In the described embodiment of the present invention the 
orientation of the locus at each point is calculated by linear 
regression of coimected edge-points in a local neighbour- 
hood (nominally 7x7). The connectivity is guaranteed by 
including only those pixels which, in {L(ij)}, share a 
conomon label. The positional-covariance matrix gives the 
orientation of the maximum variance of the data which is 
also the orientation of the least-squares fit straight line with 
two degrees of freedom. . ^ 

With reference to the definition of weighted linear regres- 
sion, with W;fc=l: 0(ij)=Qwx« is the orientation of the locus 
of edge-points without regard to the sense of the original 
edge. The ambiguity in the sense of the locus is resolved by 
borrowing from the Sobel orientation Q(r): the direction of 
the tangent to the Sobel squared magnitude iso-contour is 
always close to the direction of the tangent to the locus of 
edge-points as determined by the Canny edge-detector. 

Therefore the orientation is adjusted to minimise the 
angular separation between 0(ij) and Q(io): 

Min {<>(«J>fi(w>fai}->U^ i° {-3-3} 
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O(iJ):=O(ij>-0(ij)+l£^ii+2n7C 

where n is chosen so that the result is in [0,2n]. 
1.2 Segmentation 

line primitives are built from the remaining edges. The . 
chains of edge-pixels are broken into piece-wise straight 
sections (FIG. 9b)y and the pixels from each section 
assembled into a line structure. 

Hiere are four steps to this procedure: 

1.2.1 Vertex Analysis 

An intrinsic image {V(i j)} is constructed with labels for 
the classification of edge-pixels in {E(ij)}. According to the 
3x3 connectivity of edge-pixels in {E(io)}, one of the 
following labels arc awarded to the corresponding point in 
{V(i j)}: ISOLATED-POINT, TERMINAL, EDGE-POINT, 
T-JUNCnON,X-JUNCriON. An example of each is given 15 
in HG. 5. 

1.2.2 Comer Detection 

Comers are detected by a refinement of the EIXJE- 
POINT labels in {V(i j)}. There are two stages: 

1.2.2.1 Curve detection 20 
An intrinsic image {U(ij)} is calculated with values 

assigned to each pixel which corresponds to an EDGE- 
POINT labeUing in {V(io)}. The pixel values in {U (i, j)} 
measure the local conformity of ec^gc-pixels in {E(i j)}, in a 
7X7 neighbourhood, to a straight-line. The positional cova- 2s 
riance matrix is calculated as described in 1.1.4.3 with all 
w,=l. 

Tbe conformity, (to a straight line) is taken to be the ratio 
of major to minor eigenvalues of the positional covariance 
matrix, ie Xj©^,. Pixels previously labelled as EDGE- 30 
POINTS in {V(i j)} are then relabelled as CURVEs, if the 
conformity is less than a threshold value: 

H V(ij>=EDGE-POINT, and XJkj>OM, then V(io> 
CURVE. 

1.2.2.2 Non-maximal supression 35 
The pixels in {U(i j)} are identified as comers if they are 

both CURVE points in {V(i j)} (i.e. having a value above 
some threshold), and they are locally maximum in a 3x3 
neighbourhood. The corresponding pixels in {V(i j)} are 
relabelled as CORNERS. 40 

1.2.2 Extraction of connected edges 

Those pixels in {V(io)} still retaining an EDGE-POINT 
classification axe candidates for incorporation into the low- 
levd primidves. Such pixels are extracted from {V(ij)}. 
into {S(ij)} say, which is the linear segmentation as an 45 
intrinsic image: each connected chain of pixels is given a 
unique label in {S(iJ)}. Each cotmected set of pixels is then 
gathered into a line structure which is the low-level primi- 
tive. 

1.2.2.1 Rejection of short lines 50 

Studies have shown that sides of the signal/shadow 
boundary generally have at least some portion of their length 
for which the edge-detector furnishes an unbroken chain of 
edge-pixels. In order to avoid false positives therefore, an 
implicit rule is that a side cannot be composed entirely of 55 
very short edge-chains, but must have at least one line (in the 
line-group) which is of substantial length. Prior to segmen- 
tation and extraction of the line primitives therefore, chains 
of edge-pixels shorter than some threshold length are 
rejected (from {S(ij)}) at this point By this means, the 60 
number of low-level primitives is reduced without compro- 
mising the accuracy: short lines may form part of the final 
prospective signal-shadow boimdary following data-re- 
trieval. 

1.2.3 Weighted linear regression 65 
A weighted linear regression is performed on each of the 

lines, by which means the optimal parameters (normal 
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radius to origin, and orientation), and conformity are 
recorded. The weights are the Sobel-squared magnitudes of 
the contributing edge-pixels, as a fraction of the maximum 
Sobet-squared magnitude: 

HiCv) = G(v7/Max {G(«.v)} Vy 

1.2.3.1 Rejection of poor lines 

Lines are no longer *local* entities. The test on conformity 
to a straight line is now repeated with a less tolerant 
threshold in accordance with the extra precision afforded by 
the (arbitrarily large) number of edge-pixels in the line: 
L2.4 Linear extr^dation 

The surviving lines are extrapolated from each end, until 
they intersect either die edge of the image, or an edge-point 
which docs not belong to that line (as determined from the 
labels in {S(i j)}). The 'limits' so determined, are saved in 
the line structure. A distinction is made between each of the 
two limits which correspond to the two ends of the line. This 
distinction respects the sense of the line: see FIG. 6. 
1.2.5 Data-retrieval 

Edge-points are now recovered from {E(ij)} which were 
not classified as EDGE-POINTs in {V(ij)} (perhaps 
because diey were reclassified as CORNERS, CURVEs or 
JUNCTIONS etc). For each line, two rectangular areas are 
scanned in {E(ij)} each of which is centred on. (aligned 
wdth) the projection of the line as given by the regressed 
parameters R^lr 2nd Qw-iR- The rectangles extend from the 
line end-points, to the linuts as determined by extrapolation. 
The width of the rectangle is 5 pixels. All edge-points in 
{E(i 0)} which were not already attributed to the line are now 
included. 

Step ii) Extraction of intermediate-level primitives 
2.1 Lhiear agglomerative clustering 

In this connection the following definition are given: 
"Linear agglomerative clustering" generally is an iterative 
process for grouping togetiier lines which arc deemed (by 
some criteria) to be sufBciently similar. E.g. a pair-wise 
agglomerative clustering process based on an analysis of 
variance test, and broadly as described in "Unsupervised 
Estimation of Polynomial Approximations To Smooth Sur- 
faces in Images or Range Data** by J. F. Silverman and D. B. 
Cooper, Proa IEEE Cont on Autonaation and Robotics, 
1986, page 299^304. and defined hereinbelow, but with the 
mergable entities being lines rather than regions, and the 
constraint that two mergable regions be touching replaced 
by a constraint that two lines (or line-dusters) be lineariy 
coherent. 

"Analysis of variance" is a statistical method for com- 
paring (the similarity of) two sample distributions under the 
assumption that die underlying distribution is Gaussian; but 
the variance is unknown. The method determines an F-sta- 
tistic, which is die ratio of the two sample variances, 
correctiy normalised to take account of the true degrees of 
freedom. 

When die test is to determine the likelihood of two 
distributions having the same parent distribution, then the 
ratio is between the unexplained and the explained vari- 
ances. When applied to pair- wise clustering, the unexplained 
variance is the extra that is incurred by the proposed merge: 

where a,b arid m denote the separate, and merged entities 
respectively; % and v are the chi-squared variates and 
corresponding degrees of freedom. Hie explained variance 
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is that prior to the proposed merge: 
and 

is the F statistic, having degrees of fteedom 

V«=*V„-Va-Vil 

and 

The F statistic is used to index the ctunulative distribution 
function QOFlVo» v^) which gives the probability of that 
value or greater occuring by ctiance. If the unexplained 
variance (as a fraction of the explained) is low» thki F is 
small, and QCFtv^, v J) is large. In that case, the entities may 
be safely clustered. 

In the described embodiment of the present invention 
lines undergo linear agglomerative clustering according to a 
test of their linear coherency. For all line or line-cluster 
pairs, if the pair passes the test of linear coherency, an F 
statistic based on an analysis of variance, is computed The 
pair with the highest probability Q(Ffv^, is merged The 
parameters and limits of the merged line-cluster are recom- 
puted using weighted linear regression. Tlie whole cycle is 
repeated until no line-dusters pass the linear coherency test 
An example of a set of line-dusters, following linear 
agglomerative dustering. is given in HG. 9c. 

2.1.1 Linear Coherency 
Definitions: 

Linearly coherent l\vo lines or line-dustdrs are said to be 
lineariy coherent if they could potentially form part of the 
same side of the rectangular signal/shadow boundary. Linear 
agglomerative dustaing: an iterative process for grouping 
together lines which are deemed (by some criteria) to be 
suffidently similar. £.g. a pair-wise agglomerative cluster- 
ing process based on an analysis of variance test, and 
broadly as described in (Silverman & Coop^), but with the 
meigable entities being lines rather than regions, and the 
constraint that two'mergable regions be touching replaced 
by a constraint that two lines (or line-dusters) be linearly 
coherent 

In the described embodiment of the present invention the 
'J linear coherency lest is applied to each line or line-cluster 
pair during linear agglomerative dustering. For two line or 
line-clusters to be viable candidates for dustering: 
L They must be approximatdy colinear. 

12. Their extrapolations must overia^ or touch. 
2.1. LI Colineaiity 
Two line-clusters are determined as^Unearjg^ 
1. The Hne-clusters are approximately, paralld: 

and 

2. Th e distancc_of clo sest approach between the two 
liS^lusters is less thatTSD me - thmh old: dj.<4 pixds 
where the d^ are the distances depicted in FIG. 7. 

L2. 1.1.2 Overlapping lines 
Two line-clusters are said to have overlapjnng extre^ola- 
tions if either has a perpendicular projection onto the other, 
see nG. 8. 

2. 1 .2 Analysis Of Variance 

The F-statistic based upon the analysis of variance is 
applied to test the likdihood that two sample distributions 
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come from the same parent distribution. With reference to 
the definition of analysis of variance, and weighted linear 
regression, the dri-squared variates are the sum of square 
errors of the best fit straight line: 

where the Une-duster index j is any of a,b,m. The corre- 
sponding degrees of freedom are: 

10 

15 where and N^, arc the numbers of edge-pixels in the lines 
(of labels a and b respectively). 
Step iii) Building hypotheses 

3.1. Search Tree 

A search tree is constructed of all legal combinadons of 

20 line-dusters: see FIG. lOa-^. The tree starts with a root node 
which corresponds to the null prospective signal-shadow 
boundary that there is no coUimation material present artd 
therefore that the signal/shadow boundary is absent. The 
addition of a new branch & node to the tree represents the 

25 addition of a further line-cluster to the prospective signal- 
shadow boundary inherited from the parent node: each new 
node is assodated with a new line-cluster. Therefore the 
retum path from a particular node to the root determines the 
line-dusters in the prospective signal-shadow boundary 

^ assodated with that node. 

The search tree has at most five levels, corresponding to 
the null hypothesis, and four levels for each of the four sides 
of the rectangular signal/shadow boundary. Each node rep- 
resents a prospective signal-shadow boundary in its own 

35 right 

As each node is added, the intermediate-level tests are - 
applied to test the legality of the hypothesis. If the results are 
negative, then the node is flagged as illegal, excluded from 
further consideration, and the creation of offspring prohib- 
^ ited Assuming all the hypotheses shown are legal, the 
example of FIG. 10 would result in a search tree in its final 
state as given in FIG. 11. 

3.2. Intermediate-levd tests 

The following must be satisfied by each prospective 
*5 signal-shadow boundary during the construction of the 
search tree: 

3.2.1. Orthogonality 

All sides in a prospective signal-shadow boundary must 
be ^}proximately mutually orthc^onal. Following additicm 
50 of a new line-duster at node n, this test is applied as: 

Min{Kfi^ W-^^wijr-Jbiy2l}<0.13jk in {--4..-l.l_.4} 

where the Q are in radians, and where p ranges over all 
55 nodes but n in the return path. If the new node fails this test,^ . 
it is declared illegal. 

3.2.2. Closure 

All sides in a prospective signal-shadow boundary must 
meet, at each end, dther another side, or the edge of the 

go image. Following the addition of a new line-duster at node 
n, this test is applied by examining the limits of the new 
line-duster at node n, and other line-clusters (at nodes p) in 
the retum path which form an adjacent side of the rectangle. 
Step iv) Testing and costing hypotheses 

65 Once the prospective signal-shadow boundary has been 
built, the intersection points for each of the sides of the 
candidate signal/shadow boundary are computed Tliese 
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intersections may be either with the edge of the image, of 
with another side; see HG. 9. The high level rules must then 
be satisfied in order that the prospective signal-shadow 
boundary is verified. 

4.1. Definition of terms and symbols ^ 

Image dimensions^Let the number of pixels in the image 
array be Imagewidth and Imageheight in the x and y 
directions respectively. 

Number of sides— Let the number of (visible) sides in the 
hypothesis be N (0^Ni4). 

Intersection points — ^Let the intersection points of the kth 
side be (ajt, Pjt) and (y*. 5jt) for the intersections nearest the 
lower, and upper limits, and respectively. 

Interval — ^Lct the interval of the ktii side be 
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and 
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where tiie labels K=1.23.4 correspond to those lines 
having angles with respect to the x axis, closest to ti/2, 
Tt, 37C/2, and 2n radians. 

BOUNDARY HAS STRONG GESTALT— Tbe boundary 
has strong gestalt if: 

N=4 or BOUNDARY HAS TWO AUGNED SYMMET- 
RIC SIDES 

BOUNDARY HAS FOUR AUGNED AND ROUGHLY 
SYMMETRIC SIDES— True if: 

N = 4 



RECTANGLE NOT SMALI^Let the signal/shadow 
boundary not be small if Ajt>=10 for k in {L.N} 

Accumulated evidence — ^Let the accumulated evidence 20 
for the kth side be Wj^X^w^, j where the sum is over all 
edge-points contributing to the kth side. Let the accumulated 
evidence for the ktii interval be W^=^jW^j where the sum is 
over all edge-points contributing in the interval of the kth 
side. 25 

EVIDENCE WELL USED— Let the amount of 'unused* 
evidence be small if: WJWf^^A for all k in {I..N) 

Exposure statistics — ^Let maxjt and mean^ be the maxi- 
mimi and mean exposures (in the X-ray image) in the 
half-plane exterior to the kth side of the rectangle. Let max^ 30 
and mean^ be the maximum and mean exposures in the 
interior (signal) part of the rectangle. 

STRONG CONTRASTBOUNDARY— Let the boundary 
have strong contrast if: (maXjt-Min{X(i j)}/max^-Min{X(i, 
j)})<0.65 and meanjj/meanf, <L05 for all k in {1..N} 35 

VERY STRONG CONTRAST BOUNDARY— Let tiie 
boundary have very strong contrast if: 

(maXjt-Min{X(ij)})/(max^-Min{X(ij)})<:0.3 and 
mean^/meani„<L05 for all k in {l..>r} 

Line support— let pix^ be the number of pixels in the ^ 
interval of the kih side, for which some evidence exists, ie 
for which WjjXD. 

PERSISTENT BOUNDARY— Let the boundary be per- 
sistent if: (pix;t-H2)/Ap>0.8 for all k in 

VERY PERSISTENT BOUNDARY— Let the boundary *5 
be very persistent if: (pbcjk+2)/Ap»0.88 for all k in {l.J^ } 

BOUNDARY HAS TWO ALIGNED SYMMETRIC 
SIDES — ^Let the boundary have two aligned symmetric 
sides if: 

■ N=2 ^° 



and 



and 



(ma 
niii{-I 



.0,1} 



{!C*^-Aji/2-2nTi))s 



0.087 I >rin{1.4} 



and 

^wLR - i?^?(!zjil/Imagewidtii < 0.2 
and 

^^wLR - ^iii^l/Imagewidth < 0.21 
BOUNDARY HAS TWO LONG ASPECT SYMMET- 
RIC SIDES— True if: 
N=2 

and Imageheight ^ Imagewidth 
and 



{IQW-Kn/2-2itid} 5 0.087 ) KinU;^} 



Min 

nin {-1,0.1} 
and 

^wui - i2^{ij?l/Imagewidth < 0.1 

or 

N = l 

and Imagewidth < ImnEP^rieht 
and 



Min 

n in {-1,0,1} 
and 



{IGt«J-irjE/2-2iocl}^ 



0.087 I 



Jtin{2.4} 



{Min I 
. , in,t {lfiiai-X?i/2-2n3tl} 50,087 I . Kin {1,3} 
n in |— 1,0,1 / I 



and 

^^WLR ~ J^wLR^/l^^nBg&Nldlii < 0.1 
or 

N=t 
and 



(Min 
n in {-1.0,1} 
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^wLR - /J^^^LRl/Imagewidth < 0.1 
BOUNDARY HAS VERY STRONG GESTALT— True 

if: 

BOUNDARY HAS FOUR AUGNED AND ROUGHLY 
SYMMETRIC SIDES or BOUNDARY HAS TWO 
LONG ASPECT SYMMETRIC SIDES 
4.2. Verification rules 

A prospective signal-shadow boundary is accepted if the 
following composite nile is satisfied: 

RECTANGLE NOT SMALL and EVIDENCE WELL 
USED and (BOUNDARY HAS VERY STRONG 
CONTRAST and BOUNDARY IS PERSISTENT 

or BOUNDARY HAS STRONG CONTRAST and 
BOUNDARY IS VERY PERSISTENT 
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or BOUNDARY HAS VERY STRONG GESTALT and 

BOUNDARY IS PERSISTENT 
or BOUNDARY HAS VERY STRONG GESTALT and 

BOUNDARY HAS STRONG CONTRAST 
or BOUNDARY HAS STRONG GESTALT and 
BOUNDARY IS VERY PERSISTENTT) 
4,3. Costing hypotheses 

If a prospective signal-shadow boundary is rejected, it is 
awarded (effectively) an infinitely large cost. If the prospec- 
tive signal-shadow boundary is accepted however, it is 
awarded a cost which is determined as described below: 
Let P^ur^O.l be the fraction of pixels on the collimation/ 
signal boundary which do not have a corresponding 
edge-pixel as determined by the Canny edge-detector. 
Let P^=0.1 be the fraction of pixels which are co-linear 
with, but not on, a side of the coUimatiail/signal bound- 
ary, yet which are labelled as edge-pixels by the Canny 
edge-detector. 

Then the probability P^^ of having k or more edge-pixels 
missing on a collimation/signal boundary of total interval 
(between intersects) n pixels is : 
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where (^1-?^; Vi=2(n-4c+l); V2=^; F=4c(l-p J/(n- k-Hl). 

The same formula applies to the probability P^,^ of there 
being k or more edge-pixels which are colinear with, but not 
on, sides of the collimation/signal boundary, of total interval 
between intersects n, and total interval between limits of L. 
provided: 

The cost awarded is then: 

This function constrains the contribution from missing, 
and extra, edge-pixels to lie in [0. 0,999]. The main part of 
the cost is determined by the number of sides: 0 for 4 sides, 
1 for 3 sides, etc. The implicit rules area that: 

1. Validated hypotheses are unconditionally cheaper for 
having more sides. 

2. Given two prospective signal-shadow boundary of 45 
equal number of sides, the cheapest is that which best 
explains the data. 

Step v) Chooang the best prospective signal-shadow 
boundary 

The solution search tree is maintained by a branch and SO 
bound function. All possible hypotheses (potentially legal 
combinations of line-clusters) are considered The prospec- 
tive signal-shadow boundary system is that with the lowest 
cost, as awarded at step iv). If no other hypothesis is 
validated, the null prospective signal-shadow boundary that 
there is not signal/shadow boundary (and therefore no 
shadow region) is always vaHd (with a cost of 4 units). 

FIG. 12 shows an embodiment of an exposure and read- 
out apparatus. A radiation image of an object was recorded 
on a photostimulable phosphor screen (6) by exposing (3) 60 
said screen to X-rays transmitted through the object. The 
irradiation field emitted by the X-ray source was limited by 
using a opaque X-ray shield (4). Hie radiation image stored 
in the photostimulable phosphor screen was read-out by 
scanning the phosphor screen with stimulating rays emitted 65 
by a laser (9). The stimulating rays were deflected according 
to the main scanning direction by means of galvanometric 
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deflection device (8). The secondary scanning motion was 
performed by transporting the phosphor screen in a direction 
indicated by arrow (7) perpendicular to the main scanning 
direction. Vog stimulated emission was directed by means of 
light collector (10) onto a photomultiplier (11) where it was 
converted into an electric signal representation. The analog 
output signal is then digitised (not shown). i 

This signal can then be subjected to a processing method 
according to the present invention for the purpose of deter- 
mining the location of the signal-shadow boundary, and 
hence also to the determination of the useful image region. 
Contrast enhancing processing can then be limited to the 
diagnostically useful unage pan. 

We claim: 

1. A method for determining the signal/shadow boundary 
in an X-ray image represented by a digital signal represen- 
tation comprising the steps of: 

i) extracting lines from said x-ray image, 

ii) applying linear agglomerative clustering to said lines 
so as to form a number of line groups, 

iii) forming a number of chains of at least two of said line 
groups by appending to a single line group an addi- 
tional line group, said process of appending beiiig 
restricted by the outcome of inlermediat level tests 
performed on said single line group and said additional 
line group, 

iv) forming prospective signal/shadow boundaries, each 
prospective signal/shadow boundary comprising a 
single line group or a chain, 

v) performing high level verification tests on each pro- 
spective signal/shadow boundary and assodating a cost 
with each signal/shadow boundary, the cost being high- 
est if a prospective signal/shadow boundary comprises 
a single line group, and decreasing monotonically as a 
prospective signd/shadow boundary comprises more 
line groups, 

vi) rejecting or accepting a prospective signal/shadow 
boundary at some cost, 

vii) if at least one accepted prospective signal/shadow 
boundary exists, selecting a prospective signal/shadow 
boundary with the least cost as the signal/shadow 

' boundary; otherwise concluding that no signal/shadow 
boundary ^sts. 
Z A method according to claim 1 whereha in the linear 
agglomerative clustering process any two lines or two line 
groups, or a line and a line group are merged into a single 
line group if: 

their difference of orientation is less than some predeter- 
mined value, 

their distance of closest approach is less than some 
predetermined value, 

they have overlying extrapolations 

the outcome of the.F-statistic based upon an analysis of 
the variance of pixel distributions of lines or line 
groups to be meiged and of a prospective signal/ 
shadow boundary under consideration is such that said 
lines or line groups most likely belong to a linear pixel 
distribution. 

3. A method according to claim 1, wherein said high-level 
verification tests for tiie acceptability of a signal/shadow 
boundary, comprise a logical combination of the following 
tests for 

boundary closure, for whidi test to be passed, for a 
prospective signal/shadow boundary under, consider- 
ation, line groups are linear extrapolated at their 
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extremities so as to form extrapolated line groups and 
extrapolated line-groups must each terminate on, or 
intersect with, another line-group in the same prospec- 
tive boundary, or with the edge of the image, 

boundary shape, for which test to be passed, for the 
prospective signal/shadow boundary under consider- 
ation, is to be rectangular, or nearly rectangular, or in 
the case of there being less than four sides in a 
prospective signal/shadow boundary, the visible sides 
of said boundary are consistent with a partially visible 
rectangular or nearly rectangular shape, 

signal/shadow contrast, for which test to be passed, for the 
prospective signal/shadow boundary under consider- 
ation, is to be largely consistent with there being a 
lower level of exposure in the shadow region than the 
level of exposure in the signal region, 

boundary location, for which test to be passed, for a 
prospective signal/shadow boundary under consider- 
ation, the centre of the image must be either contained 
within the signal region, or must be close to a signal/ 
shadow boundary, 

boundary alignment, for which test to be passed, for a 
prospective signal/shadow boundary under consider- 
ation, each visible side of the rectangle must be parallel 
to a side of the image, 

boundary symmetry, for which test to be passed, for a 
prospective signal/shadow boundary under consider- 
ation, there must be two or four visible sides of a 
rectangular, which are parallel to a side of the image 
and approximately equidistant from the edge of the 
image, 

boundary completion, for which test to be passed, for a 
prospective signal/shadow boundary under consider- 
ation, the fraction of the number of pixels in an edge 
image relative to the maximnm possible number of 
pixels along the signal-shadow boundary exceeds a 
predetermined threshold. 

4. A method according to claim 1 wherein said prospec- 
tive signal/shadow boundaries comprise at least four line 
groups. 

5. A method according to claim 1 wherein said cost of 
each prospective signal/shadow boundary is further refined 
by decreasing the cost according to the fraction of the 
prospective signal/shadow boundary of which the edge 
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strength of imderlying pixels is greater than a first prede- 
termined threshold, is below a second predetermined thresh- 
old. 

6. A method of processing a digital image representation 
of an x-ray image comprising the steps of: 

(i) detennining a signal/shadow boimdary within said 
image by: 

a) extracting lines from said x-ray image, 

b) applying linear agglomerative clustering to said lines 
so as to form a number of line groups, 

c) forming a number of chains of at least two of said 
line groups by appending to a single line group an 
additional line group, said process of appending 
being restricted by the outcome of intermediate level 
tests performed on said single line group and said 
additional line group, 

d) forming prospective signal/shadow boundaries, each 
prospective signal/shadow boundary comprising a 
single line group or a chain, 

e) performing high level verification tests on each 
prospective signal/shadow boundary and associating 
a cost with each signal/shadow boundary, the cost 
being highest if a prospective signal/shadow bound- 
ary comprises a single line group, and decreasing 
monoionically as a prospective signal/shadow 
boundary comprises more line groups, 

f) rejecting or accepting a prospective signal/shadow 
boundary at some cost, 

g) if at least one accepted prospective signalVshadow 
boundary exists, selecting a prospective signaV 
shadow boundary with the least cost as the signaV 
shadow boundary; othenvise concluding that no sig- 
nal/shadow boundary exists, and 

(ii) restricting image processmg to pixels within said 
signal/shadow boundary. 

7. A method of detennining a signal/shadow boundary 
according to any of the preceding claims wherein said digital 
image representation of an x-ray image is obtained by 
scaiming a photostimulable phosphor screen that has been 
exposed to an x-my image witii stimulating irradiation, 
detecting light emitted upon stimulation and converting the 
detected light into a digital signal representation. 
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